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Recently many authors have performed sheU-model calculations of nuclear matrix elements de- 
termining the rates of the ordinary muon capture in light nuclei. These calculations have employed 
well-tested effective interactions in large scale shell- model studies. For one of the nuclei of interest, 
namely ^*Si, there exists recent experimental data which can be used to deduce the value of the 
ratio Cp /Ca by using the calculated matrix elements. Surprisingly enough, all the abovementioned 
shell-model results suggest a very small value (~ 0) for Cp/Ca, quite far from the PCAC prediction 
and recent data on muon capture in hydrogen. We show that this rather disturbing anomaly is 
solved by employing effective transition operators. This finding is also very important in studies of 
the scalar coupling of the weak charged current of leptons and hadrons. 
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The calculation of the nuclear matrix elements involved 
in the ordinary (non-radiative) capture of stopped nega- 
tive muons by atomic nuclei has been of considerable in- 
terest because they enable to access the structure of the 
effective vsreak baryonic current. Due to the large mass 
of the captured muon the process involves a large energy 
release (roughly 100 MeV) and thus surveys the baryonic 
current deeper than ordinary beta decay or electron cap- 
ture. In particular, the role of the induced pseudoscalar 
coupling Cp becomes prominent. Furthermore, the same 
nuclear matrix elements can be used in the context of 
studies of the fundamental structure of the weak charged 
current of leptons and hadrons, in particular concerning 
contributions coming from the yet undetected scalar cou- 
pling of the current. 

In the past there have been many calculations of 
nuclear matrix elements involved in the muon-capture 
processes. These calculations have been either very 
schematic ones [pj-pf or more realistic ones using various 
nuclear models ppH. Ultimately, all these calculations 
have aimed at predicting the ratio Cp/Ca of the induced 
pseudoscalar and axial-vector coupling strengths of the 
weak baryonic current by exploiting the scarce experi- 
mental information on muon-capture rates in light nuclei. 
These calculations seem to suggest wide ranges of values 
(4 < Cp/Ca < 39 for 0, 3 < Cp/Ca < 20 for ^^O 
and 28Si |,|, Cp/Ca > 13 for g), for this ratio. 
For reference the nuclear-model independent Goldberger- 
Treiman value is Cp/Ca — 6.8, which is obtained using 
the partially conserved axial current hypothesis (PCAC). 
It should be reasonable to assume that this relation be- 
tween the induced pseudoscalar and axial- vector coupling 
constants will survive in finite nuclei, although correc- 
tions may occur e.g., due to mesonic corrections in the 
weak vertices. In particular, this result should be roughly 
recovered using nuclear-structure calculations assuming 



the impulse approximation. Renormalizations within the 
impulse approximation have already been extensively dis- 
cussed for the Ca coefhcient in the context of beta de- 
cay and electron capture. On the experimental side the 
interval 6.8 < Cp/Ca < 10.6 was obtained in the or- 
dinary muon capture in hydrogen |^ and the interval 
8.8 < Cp/Ca < 10.8 was measured at TRIUMF for the 
radiative muon capture in hydrogen 

Recently, also the nuclear shell-model has been used to 
calculate the needed nuclear matrix elements for muon 
capture |@-|ll|. The calculation of was performed for 
^^C in the Os and Op shells but no definitive conclusions 
for the Cp/Ca ratio could be reached on the basis of 
their computed matrix elements. In the muon cap- 
ture transitions in ^^B and ^^B were studied but no con- 
clusions were made with respect to the possible range of 
the Cp/Ca value. In § a range of 4.1 < Cp/Ca < 8.9 
was obtained for the ordinary muon capture in ^^Na by 
fitting sevaral partial capture rates to the various final 
states in ^"^Ne. In ||ll[ several light nuclei in the Op, IsOd 
and Op-lsOc? shell-model spaces were studied by applying 
well-established effective two-body interactions in these 
spaces. In Ref. a comparison of the calculated par- 
tial muon-capture rates with the available data indicated 
a notable suppression (when compared with the above 
mentioned Goldberger-Treiman value Cp /Ca = 6.8) for 
Cp/Ca (-3.0 < Cp/Ca < 2.5) in the case of but 
nothing definitive (the interval — 7 < Cp/Ca < 15 was 
examined) could be said in case of the other studied nu- 
clei (1^0, ^°Ne, ^^Na, ^^Si and ^^S) due to controversial 
information coming from different partial capture rates 
in the same nucleus (no best fit procedure was attempted 
as in §). 

Very recently two important measurements of correla- 
tion coefficients of 7-radiation anisotropy in the capture 
of a polarized negative muon have been reported ||l^,|l2| . 
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For the allowed muon capture the angular correlation be- 
tween the emitted 7-radiation and the neutrino is scaled 
by the coefficient a [§,|o| which is related to the coeffi- 
cient 



of 11 by 



x = Mi(2)/Mi(-l) 
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Here the quantities Mi (2) and Mi(— 1) are linear combi- 
nations of the coupling coefficients of the baryonic weak 
current and the nuclear matrix elements for muon cap- 
ture [101], [121], [OUp], [nip] in the notation of 0] as 
defined in |12fl . Using the Fujii-Primakoff approximation 
(FPA) |l^] leads to the following nuclear-model indepen- 
dent result [HI 
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In the two measurements the following parameter values 
were extracted: x = 0.254 ± 0.034 [|l| and x = 0.315 ± 
0.08 0, leading to the ratios (Cp/Ca)fpa = 7.7 ± 1.2 
and (Cp/Ca)fpa — 9.7 ± 2.6, respectively, compatible 
with the PCAC prediction. 

Abandoning the FPA and determining the ratio 
Cp/Ca from Eqs. (§ and (|) using the complicated ex- 
pressions of Afi(2) and Mi(— 1), leads to rather contra- 
dictory results between different realistic nuclear mod- 
els. In [|l^ the following values of Cp/Ca were extracted 
using the different realistic nuclear models: Cp/Ca — 
3.4 ± 1.0 and Cp/Ca = 2.0 ± 1.6 for the matrix elements 
of [Q and [jl^, respectively. In the measurement of |l^] 
the corresponding extracted values are Cp / Ca — 5.3±2.0 
and Cp/Ca = 4.2 ± 2.5 [|l|. These results would in- 
dicate a small quenching of the Cp / Ca ratio with respect 
to the PCAC value. However, using the more realistic nu- 
clear matrix elements obtained from the full IsOd shell- 
model calculation with the sd-shell effective interaction 
(USD) of Wildenthal |l|, would lead to the very contra- 
dictory result of Cp/Ca = 0.0 ± 3.2 [|l^ by employing 
the computed matrix elements of Junker et a l. W% . The 
same situation occurs for the calculation of |1 l| 7verify- 
ing thereby the calculation of [|l^. This result is rather 
surprising when one takes into account the fact that the 
USD interaction is found to be very succesful for IsOd nu- 
clei in reproducing various spectroscopic quantities like 
energy spectra, Gamow- Teller decay properties, electro- 
magnetic moments and transitions (see e.g., [^|9j) as 
well as strength functions of charge- exchange reactions 

The above anomaly in the Cp/Ca predictions from 
state-of-the-art shell-model calculations is rather disturb- 
ing when contrasted with the experimental data. To give 
a deeper insight, we investigate in the present Letter the 
capture reaction f4Si(0+ ) + faAKlg 



the full shell-model framework and try to evaluate the ra- 
tio Cp / Ca through the quantity x and its measured val- 
ues reviewed above. The needed muon-capture formal- 
ism is developed in [Q[ and reviewed in the case of shell- 
model calculations in ^,|ll|. The two-body interaction 
matrix elements used in the shell-model calculation are 
given by the USD interaction [^ and a microscopic ef- 
fective interaction based on the recent charge dependent 
nucleon-nucleon interaction of Machleidt and co-workers, 
the CD-Bonn interaction [^. This is a meson-exchange 
potential model. Based on this nucleon-nucleon interac- 
tion, we derive a C-matrix appropriate for the IsOc? shell. 
This C-matrix is in turn used in a perturbative summa- 
tion of higher-order terms using the so-called Q-box ap- 
proach described in e.g., Ref. [g^. All diagrams through 
third order in perturbation theory were used to define 
the Q-box, while folded diagrams were summed to infi- 
nite order, see Ref. [|2l[[ for further details. Below we will 
label results obtained with this effective interaction by 
CD-Bonn. We have also performed similar calculations 
with the Nijmegen I [Q potential. The results were very 
close to those obtained with the CD-Bonn potential and 
hence skipped in the discussion below. 

These effective interactions are defined within the IsOd 
shell, using ^^O as closed shell core, and we perform a full 



within 



shell- model calculation using the code OXBASH 23 1. 

It is important to keep in mind that the USD interac- 
tion is fitted to reproduce several properties of IsOd shell 
nuclei, whereas the effective interaction based on the CD- 
Bonn meson-exchange potential model starts from the 
bare nucleon-nucleon interaction. The effects of the nu- 
clear medium are then introduced through various terms 
in the many-body expansion. 

The single-particle matrix elements were evaluated in 
the harmonic-oscillator basis by numerical integration of 
the radial part containing overlap of the initial and fi- 
nal harmonic-oscillator wave functions and the muonic 
s-state wave function. A harmonic oscillator basis was 
also employed in our calculation of the G-niatrix which 
enters the computation of the effective interaction, using 
an oscillator parameter of 1.72 fm. 

In addition to using the traditional bare transition op- 
erators in the evaluation of the nuclear matrix elements 
for muon capture the aim here is to employ effective 
transition operators as well. Such a calculation will be 
referred to as the renormalized one in the discussion be- 
low. To obtain effective one-body transition operators for 
muon capture, we evaluate all effective operator diagrams 
through second-order in the C-matrix obtained with the 
CD-Bonn interaction, including folded diagrams. Such 
diagrams are discussed in the reviews by Towner [ p^ and 
Ellis and Osnes |25| . Terms arising from meson-exchange 
currents have been neglected and we omit wave function 
renormalizations in the calculation of the effective op- 
erators. Intermediate state excitations in each diagram 
up to QTiLo in oscillator energy were included, as was also 
the case for the effective shell-model interaction discussed 
above. 
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As such, the determination of both the microscopic 
shell-model effective interaction and the effective one- 
body operators is done at the same level of many-body 
approach. Further details will be presented elsewhere 
pGf . In general, higher order terms introduce a correc- 
tion of the order of 10% compared with the bare transi- 
tion operator. Since the USD interaction is an effective 
one acting within the IsOd shell only, it is not possible 
to calculate a corresponding effective operator employing 
perturbative many-body methods. We will therefore em- 
ploy the effective operators obtained from the CD-Bonn 
interaction for the USD calculation as well. 

Our results are summarized in Figs. |^ and |^ with 
Cv/Ca = — 1.0 (for discussion of aspects concerning this 
ratio see for more details). 

Inspection of the nuclear matrix elements for muon 
capture of Fig. |l| shows that in the case of the USD in- 
teraction the renormalization corrections are largest for 
the largest matrix elements. This, in turn, is reflected 
in the capture rates, where the renormalized USD cal- 
culation shows much better agreement with experiment 
than the bare one. Thus, the effect of the renormaliza- 
tion is to bridge the gap between the shell-model calcu- 
lation and experiment. For the CD-Bonn interaction the 
changes in the reduced matrix elements are more subtle. 
Generally the bare matrix elements have larger absolute 
values than the renormalized ones. The USD interaction 
favours the extreme single-particle limit for the ground 
state of ^*Si, where all the valence nucleons lie on the 
d5/2 orbit. The CD-Bonn occupancies are more evenly 
distributed over the various orbits, the largest one being 
(^5/2)'^ (51/2)* (^3/2)^- Also the final state in ^^Al has 
very different structure: the USD calculation indicates 
that the partition (rf5/2)"^°(si/2)^(d3/2)"^ should be the 
most important part in the 1^ wave function whereas the 
CD-Bonn interaction predicts the state to be mostly of 
the (d5/2)'^(si/2)^(t^3/2)^ configuration. In general, CD- 
Bonn calculations favour the configurations which have 
several particles lifted up from the extreme single-particle 
shell- model limit which is, on the contrary, favoured by 
the USD interaction. 

To extract an estimate for the ratio Cp/Ca, we have 
plotted X of Eq. as a function of Cp/Ca- The ex- 
perimental value X = 0.315 ± 0.08 is taken from [ p^ . 
With this choice we obtain from Fig. || the range —0.4 < 
Cp/Ca < 1-2 for the bare USD and CD-Bonn calcula- 
tions in agreement with the USD result of ||l^ cited in 
Thus, for both of the adopted interactions the bare 
result is almost the same in spite of the basic difference in 
the origin of the used interactions and differences in the 
resulting wave functions, hinting to a strong suppression 
of the Cp /Ca ratio for the studied muon-capture tran- 
sition in the framework of the nuclear shell-model. As 
discussed before, this contradicts the shell-model calcu- 
lation of the partial capture rates in ^^Na (with a fitted 
Cp/Ca ratio) as well as experimental data on hydrogen. 

Also the USD and CD-Bonn calculations with renor- 



malized transition probabilities agree almost exactly and 
both yield a very different value for Cp/Ca than the 
calculation using bare operators. The renormalized re- 
suh is 4.4 < Cp/Ca < 5.9, close to the PCAC value 
and the results of the calculations reviewed in |l^ and 
|]lO| . This result encourages us to believe that the 
anomaly in the Cp /Ca ratio, inherent in the sophisti- 
cated shell-model calculations, has been lifted by intro- 
ducing effective renormalized transitions operators acting 
in the muon-capture process. In this way the renormal- 
ized shell- model calculations yield a value for Cp / Ca in 
^*Si compatible with data and expectations coming from 
other theoretical approaches. 

In conclusion, it is found that the renormalization 
of the transition operators is essential in the ordinary 
muon capture even in the case of an empirical effec- 
tive interaction, which otherwise reproduces the spec- 
troscopy of the involved nuclei extremely well. This is 
confirmed by the anisotropy data through the quantity 
X = A'/i(2)/Mi(— 1). The renormalization has helped 
to solve the Cp/Ca problem in the shell-model calcula- 
tions. Further analysis of this observation is in progress 
for other light nuclei in terms of the capture rates since 
anisotropy data are only available for ^*Si. For other 
nuclei of interest the theoretical situation is more com- 
plicated than is the case in ^®Si since they are either 
situated at the interface of the Op and IsOc? shells or IsOd 
and IpO/ shells, complicating thereby the evaluation of 
an effective interaction and increasing the dimensionality 
of the shell-model calculation. Research along such lines 
is in progress p^ . 
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